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Abstract 

A transition from asymmetric to symmetric patterns in time-dependent extended systems is 
described. It is found that one dimensional cellular automata, started from fully random initial 
conditions, can be forced to evolve into complex symmetrical patterns by stochastically coupling a 
proportion p of pairs of sites located at equal distance from the center of the lattice. A nontrivial 
critical value of p must be surpassed in order to obtain symmetrical patterns during the evolution. 
This strategy is able to classify the cellular automata rules -with complex behavior- between those 
that support time-dependent symmetric patterns and those which do not support such kind of 
patterns. 
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The stability analysis of patterns in extended systems has been revealed to be a difficult 
task. The many nonlinearly interacting degrees of freedom can destabilize the system by 
adding small perturbations to some of them. The impossibility to control all those degrees of 
freedom finally drives the dynamics toward a complex spatio-temporal evolution. Hence, it 
is of a great interest to develop techniques able to compel the dynamics toward a particular 
kind of structure. The application of such techniques forces the system to approach the 
stable manifold of the required pattern, and then the dynamics finally decays to that target 
pattern. 

Synchronization strategies in extended systems can be useful in order to achieve such goal. 
Different types of synchronization have been described in the literature. They have been 
applied with success to synchronize lattices of iterated maps, coupled ordinary differential 
equations and cellular automata (CA). In this note, the stochastic synchronization method 
introduced in reference for two CA is specifically used to find symmetrical patterns in the 
evolution of a single automaton. To achieve this goal a stochastic operator, below described, 
is applied to sites symmetrically located from the center of the lattice. It is shown that a 
symmetry transition take place in the spatio-temporal pattern. The transition forces the 
automaton to evolve toward complex patterns that have mirror symmetry respect to the 
central axe of the pattern. In consequence, this synchronization method can be interpreted 
as a control technique for stabilizing complex symmetrical patterns. 

Cellular automata are extended systems, discrete both in space and time. The simplest 
type of automaton is a one- dimensional string composed of N sites or cells. Each site is 
labeled by an index i = 1,...,N, with a local variable Si carrying a binary value, either 
or 1. The set of sites values at time t represents a configuration (state or pattern) a t of 
the automaton. During the automaton evolution, a new configuration <j t +i at time t + 1 is 
obtained by the application of a rule or operator $ to the present configuration: 

a t+1 = $ [a t ] . (1) 

Usually, a local coupling among the nearest neighbors is implemented. The state of site i 
at time t + 1, is a function of the value of the site itself at time t and the values of its 
neighbors s*_ x and at the same time. Thus, the evolution equation (JTJ) can be locally 
expressed as s' +1 = 0(s*_ 1; sj, sj +1 ), being the particular realization to the site level of 
the rule As there are 2 3 possible different local configurations as inputs, this means that 
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there exist 2 8 different possible evolution rules. Each one of these rules generates a different 

dynamical two-dimensional spatio-temporal pattern. Wolfram already classified the different 

Q 

structures obtained from the 256 rules in four great groups |2J. The interested reader is 
addressed to the original reference where this development can be followed. 

Our present interest resides in those CA evolving under rules capable to show asymptotic 
complex behavior (rules of class III and IV). The technique applied here is similar to the 
synchronization scheme introduced by Morelli and Zanette for two CA evolving under 
the same rule $. The strategy supposes that the two systems have a partial knowledge 
one about each the other. At each time step and after the application of the rule $, both 
systems compare their present configurations ^[(j}} and $[of] along all their extension and 
they synchronize a percentage p of the total of their different sites. The location of the 
percentage p of sites that are going to be put equal is decided at random and, for this 
reason, it is said to be an stochastic synchronization. If we call this stochastic operator T p , 
its action over the couple (^[cr*], $[of]) can be represented by the expression: 



The same strategy can be applied to a single automaton with a even number of sites. 
Now the evolution equation, a t +i = (T p o $)[at], given by the successive action of the two 
operators $ and T p , can be applied to the configuration at as follows: 

1. the deterministic operator $ for the evolution of the automaton produces $[ot], and, 

2. the stochastic operator T p , produces the result r p ($[cr t ]), in such way that, if sites 
symmetrically located from the center are different, i.e. Sj ^ sjy-i+i, then T p equals 
s N-i+i to Si with probability p. T p leaves the sites unchanged with probability 1 — p. 

A simple way to visualize the transition to a symmetric pattern can be done by splitting 
the automaton in two subsystems (cr^of), 

• a\ , composed by the set of sites s(i) with i = 1, . . . , N/2 and 

• cr t 2 , composed the set of symmetrically located sites s(N — i + 1) with i = 1, . . . , N/2, 
and displaying the evolution of the difference automaton (DA), defined as 



(2) 



(3) 
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The mean density of active sites for the difference automaton, defined as 

o N/2 
i=l 

represents the Hamming distance between the sets a 1 and a 2 . It is clear that the automaton 
will display a symmetric pattern when lim^oo p l = 0. For class III and IV rules, a symmetry 
transition controlled by the parameter p is found. The transition is characterized by the DA 
behavior: 

when p < p c —>■ lim^oop* ^ (complex non-symmetric patterns), 
when p > p c — > lim^oop* = (complex symmetric patterns). 

The critical value of the parameter p c signals the transition point. 

In Fig. ^the space-time configurations of automata evolving under rules 18 and 150 are 
shown for p < p c - The automata are composed by N = 100 sites and were iterated during 
T = 400 time steps. Left panels show the automaton evolution in time (increasing from 
top to bottom) and the right panels display the evolution of the corresponding DA. For 
p < p c , complex structures can be observed in the evolution of the DA. As p approaches its 
critical value p c , the evolution of the DA become more stumped and reminds the problem of 
structures trying to percolate the plane In Fig. E]the space-time configurations of the 

same automata are displayed for p > p c . Now, the space symmetry of the evolving patterns 
is clearly visible. 

Table D show the numerically obtained values of p c for different rules displaying complex 
behavior. It can be seen that some rules can not sustain symmetric patterns unless those 
patterns are forced to it by fully coupling the totality of the symmetric sites (p c = 1). The 
rules whose local dynamics verify (f>(si, s , s 2 ) = (fi(s 2 , s , sx) can evidently sustain symmetric 
patterns, and these structures are induced for p c < 1 by the method here explained. 

Finally, in Fig. El the asymptotic density of the DA, p t for t — * oo, for different rules is 
plotted as a function of the coupling probability p. The values of p c for the different rules 
appear clearly at the points where p — ► 0. 
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Rule 


18 


22 


30 


54 


60 


90 


105 


110 


122 


126 


146 


150 


182 


Pc 


0.25 


0.27 


1.00 


0.20 


1.00 


0.25 


0.37 


1.00 


0.27 


0.30 


0.25 


0.37 


0.25 



TABLE I: Numerically obtained values of the critical probability p c for different rules displaying 
complex behavior. Rules that can not sustain symmetric patterns need fully coupling of the 
symmetric sites, i.e. (p c = 1). 



Rule 18 Rule 150 




FIG. 1: Space-time configurations of automata with iV = 100 sites iterated during T = 400 time 
steps evolving under rules 18 and 150 for p < p c . Left panels show the automaton evolution in time 
(increasing from top to bottom) and the right panels display the evolution of the corresponding 
DA. 
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Rule 18 Rule 150 




FIG. 2: Time configurations of automata with N = 100 sites iterated during T = 400 time steps 
evolving under rules 18 and 150 using p > p c . The space symmetry of the evolving patterns is 
clearly visible. 
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FIG. 3: Asymptotic density of the DA for different rules is plotted as a function of the coupling 
probability p. Different values of p c for each rule appear clearly at the points where p — > 0. The 
automata with iV = 4000 sites were iterated during T = 500 time steps. The mean values of the 
last 100 steps were used for density calculations. 
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